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ABSTRACT 



We study a possibility to use the octopole moment of gravitationally lensed 
images as a direct measure of the third-order weak gravitational lensing effect, 
or the gravitational flexion. It turns out that there is a natural relation between 
flexion and certain combinations of octopole/higher-multipole moments which we 
call the Higher Order Lensing Image's Characteristics (HOLICs). This will allow 
one to measure directly flexion from observable octopole and higher-mult ipole 
moments of background images. We show based on simulated observations how 
the use of HOLICs can improve the accuracy and resolution of a reconstructed 
mass map, in which we assume Gaussian uncertainties in the shape measurements 
estimated using deep z'-band data of blank fields observed with Suprime-Cam on 
the Subaru telescope. 

Subject headings: cosmology: theory — dark matter — galaxies: clusters: general 
— gravitational lensing 
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1. Introduction 

It is now widely recognized that weak gravitational lensing is a unique and valuable tool 
to study the mass distribution of clusters of galaxies as well as large scale structure in the 
universe since it directly measures the projected mass distribution of the lens regardless of the 
physical state of the system and the nature of matter content (Bartelmann & Schneider 2001). 
In the usual treatment of the weak lensing analysis, the quadrupole moment of background 
galaxy images is used to quantify the image ellipticity. Then the lensing properties are 
extracted from the image ellipticities by assuming that source galaxies are randomly oriented 
in the absence of gravitational lensing. In practice we average over a local ensemble of 
image ellipticities to estimate the lensing properties. The local ensemble should contain a 
sufficient number of background galaxies to increase the signal-to-noise ratio of local shear 
measurements, whereas the region that contains the galaxies should be small enough to 
guarantee the constancy of the lensing properties over the region. The latter condition 
is necessary in the usual prescription for weak lensing because it is based on the locally 
linearized lens equation. On the other hand, the former limits the resolution of mass maps 
reconstructed via weak lensing techniques, which is of the order 1 arcmin in ground-based 
observations. 

Space-based high-resolution imaging surveys, such as the Cosmic Evolution Survey 
(Scoville et al. 2007) with the Hubble Space Telescope (HST) and the proposed Super- 
nova/Acceleration Probe (SNAP) wide weak lensing survey (Massey et al. 2004), will provide 
significant gains with a higher surface number density of well-resolved galaxies due to the 
small, stable Point Spread Function (PSF), which will enable high-resolution mapping of 
the lensing mass distribution down to an angular resolution of ~ Ofl. On the other hand, 
such a small PSF will allow us to resolve not only the elliptical component, described by the 
quadrupole moment, but also higher-order shape properties of background galaxy images, 
which could also carry some sort of information of lensing properties. It may be therefore 
interesting to see if such higher multipole moments of the shape are useful for the weak 
lensing analysis. 

There have been some attempts to generalize the weak lensing analysis to include higher 
order moments of the light distribution. Goldberg and Natarajan (2002) suggested that 
higher order effects in gravitational lensing, described by the third order derivatives of the 
lensing potential, can give rise to octopole moments of the light distribution for background 
galaxies. Goldberg and Bacon (2005) have further developed their approach and proposed 
a new inversion technique based on the Shapelets formalism (Refregier 2003; Refreger & 
Bacon 2003; Massey and Refregier 2005), and labeled this third order effect as the flexion 
of background images. Irwin and Shmokova (2006) developed a similar analysis method for 
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measuring the higher order lensing effects and applied this method to the HST Deep Field 
North. Recently Irwin, Shmokova, & Anderson 2006 reported on the detection of lensing 
signals in the UDF due to small scale structure using their "cardioid" and "displacement" 
techniques. Recently Goldberg & Leonard (2006) has extended our HOLICs approach to 
developed a method to correct HOLICs for the effect of isotropic Point Spread Function 
(PSF). 

In the present paper, following the flexion formalism by Bacon et al. (2006), we study a 
possibility to use higher multipole moments of background source images for the weak lensing 
analysis, and demonstrate via simulations how such higher order moments can improve the 
accuracy and resolution of a weak lensing mass reconstruction. 

The paper is organized as follows. After briefly summarizing the basis of weak lensing 
and the flexion formalism in section 2, we introduce higher multipole moments of galaxy im- 
ages in section 3. We define certain combinations of higher multipole moments as HOLICs 
(Higher Order Lensing Image's Characteristics) and establish an explicit relation between 
flexion and HOLICs. In section 4 we present simulations of a weak lensing mass recon- 
struction using mock observational data of image ellipticities and HOLICs. Finally some 
discussions and comments are given in section 5. 



2. Basis of Weak Lensing and Flexion 

In this section we briefly summarize general aspects of weak lensing and Bacon et al.'s 
flexion formalism. A general review of weak lensing can be found in Bartelmann & Schneider 
(2001), and we follow the notations and conventions therein. 



2.1. Local Lens Mapping 

The gravitational deflection of light ray can be described by the lens equation, 

P = 0-V1>(0), (1) 

where ip(0) is the effective lensing potential; ip is defined via the 2D Poisson equation as 
V 2 ip{6) = 2k(6), with the lensing convergence k. The convergence k = S m /S crit is the 
dimensionless surface mass density projected on the sky, which depends on the lens redshift 
z d and the source redshift z s as well as the background cosmology through the critical surface 
mass density 

c 2 D s 
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where Dd, D s , and Dd s are the angular-diameter distances from the observer to the deflector, 
from the observer to the source, and from the deflector to the source, respectively. If the 
angular size of an image is small enough to be able to neglect the change of the lensing 
potential ip(0), then we can linearize locally the lens equation ([I]) to have 5(3i = Aij(0)59j, 
where Aij is the Jacobian matrix of the lens equation, 

Aij " ( I! ) = 6ij ~ = ^ ~ ^ 6ij ~ lij ^ 
\ / ij 

where 7^ is the trace-free, symmetric shear matrix, 

H£-^M;^) (4) 

being defined with the components of gravitational shear 7 = 71 + 172. 



Q* - ' VZ K !ITZ ' (5) 



2.2. Gravitational Shear and Quadrupole Shape Moments 

In the usual treatment of weak lensing analysis, we use quadrupole moments Qij of the 
surface brightness distribution 1(0) of background images for quantifying the shape of the 
images: 

/ d 2 9 qi [I (fl^Afl.A flj 
!<POqj[I{0)\ 

where qj[I(0)] denotes the weight function used in the shape measurement, and = Qi — Qi 
is the offset vector from the image centroid. Then we define the complex ellipticity \ as 

_ On ~ Q22 + HQ12 ^ 
Qn + Q22 

The complex ellipticity x transforms under the lens mapping as 

w = X ~ 2g + g 2 X * m 
X 1 + 1^| 2 - 2Re{g X *) {) 
where g = 7/(1 — k) is the reduced shear and * denotes complex conjugate. In the weak 
lensing limit, we neglect 2nd order terms of g and x ; which yields x ~ X ~ 27. Assuming 
the random orientation of the source images, we average observed ellipticities over a sufficient 
number of images to obtain 

( X ) ^2g^ 2 7 . (8) 

The inversion equation from the shear map to the convergence map is obtained in Fourier 
space as (Kaiser & Squires 1993) 

fcr(fc) = kl "p ~? (MO). (9) 
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2.3. Spin Properties 

We define the spin for weak lensing quantities. A quantity is said to have spin-iV if it 
has the same value after rotation by 27r/N. Then, the complex shear 7, the reduced shear 
g, and the complex ellipticity x are & U spin-2 quantities. The product of spin- A and spin-5 
quantities has spin- (A + B), and the product of spin- A and spin-f?* has spin- (A — B). 



2.4. Flexion 

Flexion is introduced to be the third-order lensing effect responsible for the weakly 
skewed and arc-like appearance of lensed galaxies. The third-order lensing effect arises from 
the fact that the shear and the convergence are not constant within a source galaxy image. 
By taking higher order derivatives of the lensing potential ip(0), we can deal with higher 
order transformations of the shape quantities than the complex ellipticity. 

Flexion consists of four components of the third-order lensing tensor T>ijk = Aij t k = 
~i } ,ijk ( see Bacon, Goldberg, Rowe, Taylor 2005). The first flexion T is defined as 

jr = jr i+iJ r 2 = qqq*^ = |jr| e ^ = q Kj ( 10 ) 
and the second flexion Q is defined as 

G = Gi + i& = dddip = \Q\e 3 ^ = fry, (11) 

where d — d\ + 182 is the complex gradient operator, which transforms under rotation as a 
vector, & = de^, where is the angle of rotation. Thus T has spin-1 and Q has spin-3. The 
two complex flexion fields satisfy the following consistency relation: 

d*dG = ddT. (12) 

We then describe the transformation of the shape of a background source by expanding 
the lens equation (pQ) to the second order as 

d/3i « Ajddj + -V^dOjdOk. (13) 

The third-order lensing tensor T>ij k can be expressed as the sum of the two terms, T>^ k = 
Fijk + Qijki with the spin-1 part and the spin-3 part Gijk'- 

c _ 1 (Qi Q2 \ c _ 1 ( Q2 -Qi\ n ,. 
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Flexion has a dimension of length -1 (or angle -1 ). This means that the effect by flexion 
depends on the source size. The shape quantities affected by the first flexion T alone have 
spin-1 properties, while those affected by the second flexion Q alone have spin-3 properties. 

From equations ({TO]) and (iTTj) . the inversion equations from flexion to the convergence 
can be obtained as follows (Bacon et al. 2006): 

(k + iB)^ = A^JF, (16) 
(n + iB)g = A- 2 3*d*d*g, (17) 

where the complex part %B describes the B-mode component that can be used to test the 
noise properties of weak lensing data. An explicit representation for the inversion equations 
is obtained in Fourier space as follows: 

M*0 = ~% , 2 , 2 , (18) 
. .. &(fc?-3M:a)+g2(3fc 1 2 fc2-fcS) nq , 

Kg{k) ~ 1 W+W 2 ' 

for k ^ 0. Further we can combine independent mass reconstructions k a (k) (a = ^f,J-,Q) 
linearly in Fourier space to improve the statistical significance of the k map with minimum 
noise variance weighting: 

*(*) = E "^° (fc) " a(fc) , (20) 

where W Ka (k) = 1/P^\k) with noise power spectrum P£p(k) of a k map reconstructed 
using ath observable: 



Pi N \k) a 2 



7 



2 8nn g 



2k 2 %Tin g k 2 

with Pa (k) being the shot noise power of ath observable, a a being the intrinsic dispersion of 
ath observable, and n g being the surface number density of background galaxies. Assuming 
that errors in k a (k) between different observables are independent, the noise power spectrum 
for the estimator (120]) is obtained as 

Pf'(fe) = 1 = Kw^- (22) 

E a w a (k) E a yp^(k) 
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3. Higher multipole moments of images: HOLICs 

In this section we consider higher multipole moments of images and define useful combi- 
nations of them as Higher Order Lensing Image's Characteristics (HOLICs). We then derive 
a simple, explicit relation between flexion and HOLICs. 

Higher order moments of images are defined as a straightforward extension of the 
quadrupole moment. The octopole moment Qijk and the 16-pole moment Qijki are define as 
follows: 

fd?e qi [i(0)] ' 1 ] 

^ ijkl ~ fd?e qi [i(o)] ' 1 } 

We first define the normalization factor £ as 

£ = Qua + 2Q1122 + Q2222 (25) 

with spin-0. Then, we define the following combinations of octopole moments as our HOLICs: 

^ = (Qui + Q122) + i{Qin + ^ 222 ) (26) 
, {Qui - 3Q122) + «(3Qn2 - Q222) /ov x 

* = 1 ' (27) 

where the first HOLICs ( has spin-1, and the second HOLICs 5 has spin-3. Note that 
HOLICs have the dimension of [length] -1 (or angle -1 ), the same as flexion does. 

Now we are in a position to derive the transformation law of HOLICs under gravitational 
lensing. For this purpose we first derive the relation between the source octopole moment 
Q\*l and the image octopole moment Qijk- A straightforward calculation leads to 

1 

2 l 

+vAjivAfc n 'Pj mo AAn*A.j m *A.knFo)Q 'imno (28) 



Qijk ~ •Au*A.j m tA.k n Ql mn ~\~ {■A-il^-jmP kno ^-jm^kn^ilo 



where we have used the fact that the integration measures in the source and image planes 
are related in the following way (see Appendix A for detailed calculations): 



d 2 {3 



d/3 



d 2 e 



06 

(1 -k) 2 (1 -2FiA0i -2F 2 A6 2 )d 2 6, (29) 
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to the first order of reduced flexion defined as 

F := F 1 + iF 2 = -^— (30) 

1 — K 

G := G 1 + iG 2 = T ^—. (31) 

1 — K 

Note that the flexion term from the determinant does not yield a net contribution to the 
denominators of equations (|23|) and (|24|) since the coordinate system is taken such that the 
first moment of 1(0) vanishes: 

J d 2 p qi [I(6)] « (1 - Kf J d 2 9 qi [I{0)l (32) 

where we have neglected the second order term in 7. From this transformation law one 
obtains the desired expressions as 

&) C - 2gC - 9*S - z(8F* V + 9F + IGrf + G*\) 

^ (1 - «)(1 - ARe[g*X\ - 5Re[Fi}] - Re[Gt* ni \) ' 1 ' 

, W 5-3gC-i(10FT, + 7F*\ + 3G) 

(1 - k)(1 - ARe[g*X] - 5Re[Fi*j] - Re[Gt} n \) ' 1 ' 

where dimensionless quantities 77 and A are defined with 16-pole moments by 

_ (Qnu — Q2222) + 2i(Qni2 + Q1222) /„.v 

. (Qini - 6Q1122 + Q2222) + 4z((5ni2 - Q1222) ,„„,. 
A = , (36) 

with spin-2 and spin-4, respectively; t/, t/jj and iy are defined with 32-pole moments by 

(Q1IIII+2Q11122+Q12222) + ^(Qllll2+2<5ll222+Q22222) 



tv = 



(Qlllll~2QlH22~3Ql2222) + ^(3Qllll2+2Q 112 22~Q22222) 

(Qiini - IOQ11122+5Q12222) + *(5<5nii2— IOQ11222+Q22222' 



(37) 



with spin-1, spin-3, and spin-5, respectively. 



Assuming that the quantities g, rj and A are small and neglecting higher-order terms 
containing 32-pole moments which are reasonable assumptions on the weak lensing data, we 
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can approximate the above equations as 

C« « C-^, (38) 
6 (s) w (39) 

41 - AC 

The formulae (138!) make it possible to relate directly the flexion fields and the HOLICs 
measurements. Since the ( and 5 are quantities with non-zero spin, namely quantities with 
directional dependence, the expectation values of intrinsic ( and 5 are assumed to vanish, 

<o » li^"!'- < 40 > 

Neglecting the flexion term in the Jacobian matrix (|29|) will lead to a reduction of the 
response (O/J 7 from 9/4 to 5/4 while it will keep the response (8)/Q unchanged, which was 
found earlier by Irwin & Shmakova 2006. In this way, one can measure directly the flexion 
fields T{6) and G(0) from the observable HOLICs. Once we obtain the flexion fields, we can 
make use of equations (115)1 and (119)) to invert them to the surface mass distribution. 



It is important to note the above relation (1401 is modified if we take into account the fact 
that the "apparent" center of an image defined by the first moment of the image is different 
from the "actual" center mapped by the lens equation from the center of the source. We 
discuss in detail this shift of the centroid in Appendix B. 



4. Simulated Observations 

In order to test the performance of mass reconstructions based on HOLICs measure- 
ments, we generate simulated observations of the weak lensing effects, namely (x, C that 
include observational errors as Gaussian uncertainties. The flexion fields (JF(0), Q(6)) can 
be used to reconstruct mass maps directly, independent of information on the shear field 
7(0) (see §2.4). The equation (135)) defines the direct, unbiased estimators for the flexion 
fields, where the precision of this measurement depends on the intrinsic values of HOLICs 
convolved with the measurement noise. In the present simulation we do not take into account 
explicitly the centroid shift (see Appendix B) but use directly equations ( 133)) and (l34j) to 
calculate the lensed HOLICs from the intrinsic shape quantities and lens properties, which 
does not require the measurement and removal of the apparent centroid of an image. 

To determine the widths of random Gaussian distributions for the noise component 
of HOLICs {(,6), arising from the intrinsic scatter in unlensed HOLICs and observational 
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noise, we refer to variances of HOLICs obtained from our preliminary study of deep ?'-band 
data of ~ 2 x 4 deg 2 blank fields observed with Suprime-Cam on the Subaru telescope 
(T. Yamada, private communication). Each 4deg 2 data set consists of 18 Suprime-Cam 
pointings with C/202 pixel -1 . We used our weak lensing analysis pipeline based on IMC AT 
(Kaiser, Squires, & Broadhurst 1995) extended to include the higher multipole moments 
in the shape measurements. We selected a sample of 493934 background galaxies with 
20 < i' AB < 24.5 mag in the blank fields, corresponding to a mean surface number density 
of n g ~ 17arcmin -2 . Here we discarded stars and all objects for which reliable shape 
measurements cannot be found. In particular, we excluded those small objects whose half- 
light radius (r h ) is smaller than 0'/50 (2.5 pixels). We note the median value of stellar r h 
over the entire field is (r£) mcd 0'/40 with a dispersion of 0'.'08 using iV* = 27958 stars 
(n* ~ larcmin -2 ). This lower cut-off in the galaxy size is essential for us to be able to 
make reliable shape measurements because (1) the smaller the object, the noisier its shape 
measurement due to pixelization (i.e. discretization) noise, in particular for the case of 
higher-order shape moments, and (2) the shape of an image whose intrinsic size is smaller 
than or comparable to the size of PSF can be highly distorted and smeared. For example, if 
the spatial distributions of both the source and the PSF are described by a two-dimensional 
Gaussian, with half-light radii of r^o and r£, respectively, then the half-light radius of the 
PSF-convolved image is r h = \Jr\ Q + r* h 2 . When both the sizes are equal, then r h = V2r* h , 
which is close to our choice for the lower cut-off in rv 

We estimated the unlensed dispersions of HOLICs to be = a/ (|C| 2 ) ~ 0.0215 arcsec -1 
and a 5 = y/(\5\ 2 ) ~ 0.0248 arcsec , where noisy outliers responsible for the non-Gaussian 
tail were removed in the variance estimation. After clipping rejections the number of galaxies 
in this "clean sample" is about N c i ean m 3.1 x 10 5 (n g ~ llarcmin -2 ), and the median value 
of rh is (r£) ~ 0'/639. Using the clean sample we also measured dispersions of dimensionless 
HOLICs, a( and aS, with a being the characteristic scale of the observed galaxy (Goldberg 
& Bacon 2005; Goldberg & Leonard 2006). We take a to be the half-light diameter, a = 2r^, 
while Goldberg & Bacon (2005) and Goldberg & Leonard (2006) chose a to be the semimajor 
axis for measuring the intrinsic flexion of galaxies. We found a a ^ = 0.0286 and a a & = 0.0344 
with the median size of (a) me d = 1"27. We note that there is a good agreement between 
°af/ ( a )mcd (o'as/ (a)mcd) and o"^ (erg), indicating that the measured dispersions are effectively 
weighted by moderately large galaxies with a xs V . The first HOLICs (, which is a spin-1 
quantity, is sensitive to the determination of the centroid, or the first moment, with spin- 
1. Thus the error in the centroid determination affects seriously the estimation of the first 
HOLICs (, while this effect is of second order for other shape quantities with non spin- 
1 quantities. Note that even though small objects were excluded from the analysis, no 
anisotropic/isotropic PSF corrections were applied in measuring the HOLICs for the present 
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study, which implies that the dispersions of HOLICs were underestimated to some degree 
due to the isotropic smearing effect by the atmosphere and the Gaussian weighting in q[I(0)] 
(see Fig. 2 of Irwin et al. 2006). 

In the present study, we simply adopt for illustration purposes random Gaussian dis- 

(s) (s) 1 

tributions with dispersions of = <j s = 0.02 arcsec in generating the intrinsic val- 

(s) i (s) 

ues of HOLICs, corresponding to flexion dispersions of a F ~ 0.009 arcsec and a G ~ 
0.027 arcsec -1 . We note that Bacon et al. (2006) adopted = (Jq* = 0.04 arcsec -1 in 
their simulation. For the intrinsic dispersion of ellipticities, we adopt = 0.4. As a lens 
model we use a mock cluster located at z& = 0.1, which is shown as the dimensionless surface 
mass density k in the top- left panel of Fig. 1. The cluster consists of a main halo and five 
sub-halos which are described by NFW density profiles with different masses and concen- 
tration parameters. The field is square-shaped with a side length of 6.4 arcmin. We use 
N s = 4, 096 sources {n g = 100 arcmin -2 ) randomly distributed over the field, with intrinsic 
shape quantities (x , C j drawn from random Gaussian distributions with dispersion 
(o"x \ a [ S \ a )> respectively. We assume the background sources are located at a single red- 
shift of z s = 1.2, which is a fair approximation for the lens located at a low redshift of 
Zd = 0.1. These observation parameters are appropriate for a future space-based survey such 
as the planned weak lensing survey with the SNAP satellite (Massay et al. 2004), except 
that the intrinsic flexion dispersions (<Jp\ Oq) are derived from ground-based weak lensing 
data with galaxies of a = 2r^ m 1". Hence the values of (<7p \ <j^q) adopted in the present 
study will probably be very optimistic for space-based data with distant galaxies of a ~ 0f!5 
(see Fig. 2 of Massey et al. 2004)0 Such small, distant galaxies will have intrinsic flexion 
dispersions greater by a factor of a few. Finally, we use equations d7j), (I3"3"j) . and (1541) to 
generate lensed quantities (x, C, <5) for all sources. 

Figure 1 shows reconstructed mass maps of the model cluster using (x, C, S) along with 
the input mass model. We used the linear inversion equations (Q, ffl8l) . and ffl9|) . The recon- 
structed /t-maps were smoothed with a Gaussian filter, where the Gaussian FWHMs are taken 
to be (0'333, 0!083, 0'249) for reconstructions using (7, F, G), respectively. Then, the disper- 
sions in the Gaussian-smoothed 5-mode maps are obtained as ((Tb(j), (Tb(F), (t_b(G)) = 
(0.0284,0.0236,0.0570). The k maps reconstructed using HOLICs recover substructures 
better than large-scale structures, allowing a high-resolution mass reconstruction. The su- 
perior sensitivity of flexion to small scale structure comes from the l//c 2 -dependence of the 
noise power spectrum in a k reconstruction (see equation [21]) . Further, since we have as- 
sumed ai = 0£ whereas ( has a three-times larger response to flexion than 5, the JF-based 



1 For a galaxy with a Gaussian profile, a = 2r/j = VS In 2a « 2.354cr with a being the Gaussian dispersion. 
Hence, a ~ 0'.'5 for a - Of! 2. 
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reconstruction has a three-times better sensitivity than the (/-based reconstruction. 

In Fig. 2 we show the k map obtained by combining shear- and flexion-based recon- 
structions using equation (|2"U1) . On small angular scales the signal is dominated by flexion, 
while it is dominated by gravitational shear on large angular scales. Since W K p{k) > W K a{k) 
for all wavenumbers, the k map is less weighted by the second flexion for the observation 
parameters adopted in this study (i.e., = a§). The rms noise level in the K-map from 
flexion and shear data is reduced down to 

a B « 1 = 0.0173. (42) 

5. Discussion and Conclusions 

In the present paper, we have studied the possibility to improve the weak lensing analysis 
by utilizing the octopole and higher-multipole moments of lensed images that carries the 
third-order weak lensing effect. By defining proper combinations of octopole moments as 
HOLICs, we have derived explicit relations between the flexion fields (F, G) and observable 
HOLICs (£, 5). In the weak lensing limit, the first flexion F excites in lensed images the first 
HOLICs ( with spin-1, while the second flexion G excites the second HOLICs 5 with spin-3. 
One can employ the assumption of random orientation for intrinsic HOLICs of background 
sources to obtain an unbiased, direct estimator for flexion, in a similar manner to the usual 
prescription for weak lensing. 

We have also shown by using simulated observations how the use of HOLICs can improve 
the accuracy and resolution of a reconstructed mass map, in which we assumed Gaussian 
uncertainties in the shape measurements estimated using deep i'-band data of blank fields 
observed with Subaru/Suprime-Cam. The gravitational shear and flexion have different 
scale-dependence in mass reconstruction errors. The mass maps reconstructed using HOLICs 
recover substructures better than the shear-based reconstruction, allowing a high-resolution 
mass reconstruction. It is shown that an optimal linear combination of individual mass 
reconstructions can be formed using the statistical weight in Fourier space, which can improve 
the statistical significance of weak lensing mass reconstructions. In actual observations, on 
the other hand, we must apply various shape corrections (e.g., isotropic/anisotropic PSF 
corrections) to the higher order shape quantities in order to measure flexion to high precision. 
In particular, the first HOLICs £ with spin-1 is highly sensitive to the choice of the center- 
of-image. These issues will be discussed further in the forthcoming publications. 
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Appendix A: Calculation of the Jacobian 

We present detailed calculations of the Jacobian (29) in this appendix. Using the lens 
equation the Jacobian can be calculated as follows: 

(Au + V in d9i + V n2 de 2 ) (A22 + T>i 22 d6i + V 222 d6 2 ) - (Ai 2 + Vu 2 d6i + Di 22 *) 2 j d 2 9 

{l-K) 2 X 

l-g- l(3Fi + Gi)d6i - ]-{F 2 + G 2 )d6 2 \ (l + g- ^(Fi - G 1 )d9 1 - \(?>F 2 - G 2 )d9 2 



2 V ^ ^ 2 y ' 7 V 2 V 2 



1 - k) 2 (1 - 2Fid0i - 2F 2 d0 2 ) d 0. 



Appendix B: Effect of the Centroid Shift 

In a weak lensing analysis we quantify the shape of an image by measuring various 
moments of the surface brightness distribution 1(0), in which the moments are calculated 
with respect to the centroid of the image, or the center of light, defined by the first moments 
of 1(0): 

= jd 2 eej(0) 

l ~ Jd^9i(0) ■ { > 

However, in general, this apparent center can be different from the "point" that is mapped 
using the lens equation from the center of unlensed light. We refer to this point as the 
"true" center of the image. The difference between these two centers, namely the apparent 
and the true centers, cause a significant effect in evaluating the first flexion as pointed out 
by Goldberg & Bacon (2005). For the second flexion the effect of the centroid shift is second 
order so that we will ignore it in the present paper. 

Let us define the center of the source in the absence of gravitational lensing as $ and 
the apparent center of the lensed image as We also define the true center of the image 
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that is mapped from /3j using the lens equation, as 9i(/3). Similarly, we define the point in 
the source plane that is inversely mapped by the lens equation from 9^ as /3i(0). Then, we 
have 

A(0) - Pi = Oi ~ at(p) - A = 9 t - at(p) - ! f^f s ^ 

e ( n\ $d 2 9(l-2F 3 d9 j ) (9 l -a l (6))m ^ 
~ 9i - a W Jd 2 9(l-2F j d9 J )I(6) (3) 

where we have used the expression (pQ) for the Jacobian and used the relation between the 
lensed and unlensed surface brightness distributions, 1(6) = I s '((3), with $ = 9^ — «j (#,,). 
We then expand the deflection angle ai in equation (jSJ) as 

an(0) « ati{d) + a id d9j + -anjkdejdBk = ajfi) -(A- 1)^ d9 3 - - ]-V ijk d9 3 d9 k (4) 

with 9i = $i + d9{. Thus, in the first order of the gravitational shear and flexion, we have 

A(0)-A = Bi-onifi) 

J d 2 9 (1 - 2F l d9 l ) (9i + d9 t - oi(0) + (A- d9 3 + \V l]k d9 3 d9 k ^ 1(6) 

Jd 2 9(l-2F l d9 l )I(6) 
/ d 2 9 (1 - K,)2F l d9 l d9 i I(0) 1 / d 2 9V ijk d9 j d9 k I(6) 
fd 2 9I(6) 2 Jd 2 9I(0) 

= t TQf^ i + ^*x)i + l[Sx%\ (5) 

where trQ = Qn + Q22 is the trace of Q i3 defined by equation Thus the displacement 
from the true to the apparent center is trQ (|jFj + |[.F*x]i + |[£?X*]i) i n the "source plane", 
and all of the shape quantities in equation (jSJ) are measured using the apparent center 6. In 
the zeroth order of the shear and flexion, this displacement is magnified by 1/(1 — k) in the 
image plane, so that we have 

0,(0) « Si - HQ Qf 4 + \[F* X \ t + \[Gx m ]l\ (6) 

where F and G are the first and the second flexion, F = T j(\ — k) and G/(l — k). 

Next, we evaluate the effect of the centroid shift for measuring the first HOLICs, (. We 
define the apparent octopole moments (with respect to the apparent center 6) as 

Jd 2 9d9 i d9 j d9 k I(0) 
Qljk ~ Jd 2 9I(0) (7) 
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and the "true" octopole moments (with respect to the true center 6(0)) as 

/ d 2 6 (6i - Biifl)) (e 3 - 6,0)) {Ok - 0*03)) 1(6) 



^ ijk Jd 2 6I(0) 

Jd 2 e(d9 i + At)(dB j + A e j )(d9 k + Al)I(0) 



fd 2 9I(6) 

where A^ = tiQ + §[-F*x]i + |[^X*]i)- Using the expressions for the moments above, it 
is straightforward to calculate the corresponding HOLICs using the moments with respect 
to the true center, and then to have the relation between £' constructed using the true center 
and ( constructed using the apparent center: 

C = C + 2^ A* + ^ixA**} « C + ^ (3F + AF* X + l -G X ^ (10) 



as well as the relation between 5 and 5: 



S^S + S^lxA^^S + l^xF, (11) 

up to the the second order of the shear, flexion and non 0-spin shape quantities. Equations 
(flUj) and (jTTI) allow us to establish the required relationships between the unlensed HOLICs 
(C j<^) an d the lensed HOLICs (C,S) (note that equations [33] and [31] are the relations 
between and £*, and between 5^ and 5*, respectively, in the notation of this appendix), 

C - 2g(* - g*5 - \(8F*r] - 1Q^^F* X + 9F- 12^f^F + 2Gr]* - 2 { ^f^G X * + g*A) 
^ ~ (1 - k)(1 - 4Re[g*X] - 5Re[F$\ - Re[Gt} n }) ^ 

5 - 3g( - \(lOFr] - 18^^F X + 7F*X + 3G) 

(1 - k)(1 - 4Re[g*\] - 5Re[Ft}} - Re[Gt* m ]) ' V ; 

Finally, in the first order, the first flexion is expressed using the observable shape quantities 

as 

(0-(!-s(Sf!))* (14) 

One can see the asymptotic behavior of the correction term in equation (fl4l) . when the lensed 
image is close to a circle. For an image with the circular top-hat brightness profile, we have 
(trQ/£) 2 = 3/4, so that one will measure ( = regardless of the value of F. 
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Fig. 1. — Reconstructions of the dimensionless surface mass density n of a model cluster 
at Zd — 0.1 from simulated ellipticity and HOLICs data. The side length is 6.4 arcmin, and 
4,096 (n g = 100arcmin~2 ) randomly generated galaxies at a redshift of z s = 1.2 have been 
used for reconstructions (see §4). Shown in the top-left panel is the input k map of the 
model cluster. The top-right panel shows a k map reconstructed from image ellipticities, x- 
The bottom panels show k maps reconstructed using HOLICs, (c) ( and (d) S. The lowest 
contour and the contour interval are Ak = 0.06. The shaded circle in each panel indicates 
the FWHM of the Gaussian filter. 




Fig. 2. — Dimensionless surface mass density k obtained by combining shear and flexion 
data in Fourier space (see equation [20]). The lowest contour and the contour interval are 
Ak = 0.06, 



